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ABSTRACT 

We present results of the optical spectral and photometric observations of the 
nucleus of Markarian 6 made with the 2.6-m Shajn telescope at the Crimean Astro- 
physical Observatory. The continuum and emission Balmer line intensities varied more 
than by a factor of two during 1992-2008. The lag between the continuum and H/3 
emission line flux variations is 21.1 ± 1.9 days. For the Ha line the lag is about 27 days 
but its uncertainty is much larger. We use Monte-Carlo simulation of the random time 
series to check the effect of our data sampling on the lag uncertainties and we compare 
ou r simulation results w ith those obtained by random subset selection (RSS) method 
of iPeterson et al.1 (1998). The lag in the high- velocity wings are shorter than in the 
line core in accordance with the virial motions. However, the lag is slightly larger in 
the blue wing than in the red wing. This is a signature of the infall gas motion. Prob- 
ably the BLR kinematic in the Mrk 6 nucleus is a combination of the Keplerian and 
infall motions. The velocity-delay dependence is similar for individual observational 
seasons. The measurements of the H/3 line width in combination with the reverbera- 
tion lag permits us to determine the black hole mass, Mbh — (1.8 ± 0.2) x 10 8 Mq. 
This result is consistent with the AGN scaling relationships between the BLR radius 
and the optical continuum luminosity (Rblr °c L ' 5 ) as well as with the black-hole 
mass-luminosity relationship (Mbh — L) under the Eddington luminosity ratio for 
Mrk 6 to be L bol /L Edd ~ 0.01. 



Key words: 

ual: Mrk 6 



galaxies: active - galaxies: nuclei - galaxies: Seyfert - galaxies: individ- 



1 INTRODUCTION 

Over the past n early 30 years t he method of reverberation 
mapping (RM) (|Petersonlll988l and reference therein) has 
become one of the standard methods for studying the Ac- 
tive Galactic Nuclei (AGNs). It is based on the assumption 
that in a typical Seyfert galaxy the source of continuous ra- 
diation near a black hole, named as accretion disc (AD), 
is expected to be of order 10 13 — 10 14 cm. Photoionization 
of the gas located at a distance of order 10 16 cm produces 
broad emission lines. The relationship between the contin- 
uum and emission line fluxes can be represented by the equa- 
tion (JBlandford fc McKedll982h 
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L(t) = / ^(r)C(t - r)dt, 

where C(t) and L(t) are the observed continuum and 
emission-line light curves, and 9(t) is the 1-d transfer func- 
tion (TF). The TF determines the emission line response to 
a 5-function continuum pulse as seen by a distant observer. 
So, the emission lines "echo" or "reverberate" in response 
to the continuum changes with a delay r. The size of the 
region where broad lines (BLR) are formed can be written 
as Rblr = ct. The primary task of the RM method is to 
use the observable C(t) and L(t) to solve the above integral 
equation for the TF in order to obtain information about 
the geometry and physical conditions in the BLR. Unfor- 
tunately, it is very difficult to find a unique and reliable 
solution to this equation. However, it is possible to find a 
temporal shift (lag) between the continuum and emission 
line light curves using the cross-correlation analysis. 
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Applying the virial assumption, the mass of the black 
hole can be determined when the BLR size and the veloc- 
ity d ispersion of the BLR gas are known (|Peterson et al.l 
|200J). The present tremendous progress being made in black 
hole mass estimates can be attributed to the reverberation 
method. 

On the other hand, different segments of a single emis- 
sion line seem to be formed at different effective distances 
from the ionizing source. In that case, the response in the 
flux of emission line at line-of-sight velocity V r and time de- 
lay t is caused by the 2-d transfer function or "velocity delay 
map" iJHorne et al] 12004 ). The reverberation technique ap- 
plied to different parts of a single emission line allows one to 
make conclusions about the velocity field of the BLR gas. To 
the present time, considerable progress was made in under- 
standing the direct ion of the BLR gas motion (Gaskcll 2009; 
iBentz et al] |2008| ; iDennev et all 12009) : iBentz et al] l2010h . 
For example, some distinctive signatu res for infall i ng gas 
motions in NGC 3 516 and Arp 151 JBentz et al] l2009bl ; 
IDennev et al] 120091 ) were revealed: the blue side of the line 
lagging the red side. NGC 5548 shows the virialized gas mo- 
tions with the symmetric lags on both the red and blue sides 
of line. However, the BLR gas in NGC 3227 shows the sig- 
nature of radial outflow: shorter lags fo r the blue-shifted ga s 
and longer lags for the red-shifted gas (|Dennev et al.ll2009l) . 

More than 40 AGNs have been studied by the RM up to 
now (|Peterson et al.ll20oi : IBentz et af1l2009bl ; IDennev et ail 
[2010j). However, the Mrk 6 nucleus is absent in this list. 
Just a few studies of this galaxy have been made by 
the reverberation method. W e can mention the paper by 
iDoroshenko fc Sergeevl (J2003J) based on the archive spec- 
tra of Mrk 6 obtained from 1970-1991 using the image 
tube spectrograph at the 2.6-m telescope of the Crimean 
Astrophys ical Observat o ry (C rAO). There is also another 
paper by ISergeev et al] (119991 ) that includes the results of 
the 1992-1997 observations with the same spectrograph but 
with a CCD detector. In this paper a lag in the flux vari- 
ations of hydrogen lines with respect to the adjacent con- 
tinuum flux variations was reported for the first time and 
changes in the line profiles were studied. 

Mrk 6 is a Seyfert 1.5 galaxy (Sy 1.5). This galaxy 
was one of the first galaxies in which the strong vari- 
ability of the H/3 emis si on li ne profile was detected by 
iKhachikian fc Weedmanl (|197ll ). Although there have not 
been many optical observations of M rk 6, there are some 
radio studies of Mrk 6 available (IKukula et al] 1 19961 : 
iKharb et al] 120061) ■ which revealed the complex struc- 
ture of the radi o emission. The X-ray emission from the 
Mrk 6 nucleus (iFeldmeier et al.l ll999l: iMalizia et ail 120031 : 
llmmler et alj|2003l ; ISchurch et al.l l2006) exhibits a complex 
X-ray absorption and some authors assume th at the BLR is 
a pos sible location of this absorption complex (|Malizia et al] 
l2003h . 

In this paper we present the results of our optical spec- 
troscopic observations of Mrk 6 during the monitoring cam- 
paign from 1998 to 2008 performed after publishing our 
earlier papers with observations made from 1970-1997. For 
completeness we apply our analysis to all the spectral CCD 
observations that have been made since 1992 and we also 
use the results of our photometric observations in the V 
band. Throughout the entire work, we take z = 0.01865 
from our spectral estimates, the distance to Mrk 6 equal to 



D = 81 Mpc, and H =70 kms" 1 Mpc" 1 . The observations 
and data reduction are described in Section 2. We present 
the cross-correlation analysis in Section [3] the line width 
measurements in Section [4] the estimates of the black hole 
mass, the mass-luminosity and lag-luminosity diagrams are 
presented in Sections [5] and [6] the velocity-resolved rever- 
beration lag analysis is performed in Section [7] The results 
are summarized in Section [8] 



2 OBSERVATIONS 

2.1 Spectral observations and data processing 

The Ha and H/3 spectra of the Seyfert galaxy Mrk 6 were ob- 
tained from the CrAO 2.6-m Shajn telescope. Prior to 2005 
we used the Astro-550 CCD which had a size 580x520 pix- 
els and was cooled by liquid nitrogen. The dispersion was 

2.2 A pixel -1 , and the spectral resolution was about 7-8 A. 
The working wavelength range was about 1200 A. The en- 
trance slit width was 3". For technical reasons it was not 
always possible to set the same position angle (PA) of the 
entrance slit. Most of the observations were performed along 
PA~ 125° or PA~ 90°. The typical exposure time was about 
one hour. The "extraction window" was equal to 11". 

In July 2005 the Astro-550 CCD was replaced with the 
SPEC-10 1340 x 100 pixel CCD, thermoelectrically cooled 
up to — 100°C. In this case the dispersion was 1.8 A pixel -1 . 
A 3'.'0 slit with a 90° position angle was utilized for these ob- 
servations. The higher quantum efficiency (95% max.) and 
the lower read noise of this CCD allowed us to obtain higher 
quality spectra under shorter exposure times. The spectral 
wavelength range for these data sets was about 2000 A near 
the Ha and H/3 regions. However, the red and blue edges of 
the CCD frame are unusable because of vignetting. Gener- 
ally, all spectra in both spectral regions were obtained with 
a single exposure during the whole night. Our final data set 
for 1998-2008 consists of 135 spectra in the H/3 region and 
48 spectra in Ha. The mean signal-to-noise ratio (S/N) is 
equal to 32 for the Astro-550 and 73 for the SPEC-10 in 
the H/3 region and is equal to 40 and 66 in the Ha region, 
respectively. 

As a rule, each observation was preceded by four short- 
time (~10 s) exposures of the standard star BS 3082, whose 
spectra were obtained at approximately the same zenith dis- 
tance as that of the galaxy. The s pectral energy distr i butio n 
for these spectra was taken from iKharitonov et al] (|l988f) . 
The standard star spectra were used to remove telluric ab- 
sorption features from the spectra of Mrk 6, to provide the 
relative flux calibration, and to measure the seeing param- 
eter defined as the FWHM of the cross-dispersion profile 
on the CCD image. The description of the primary data 
processing as well as some details of absolu te calibration 
and m easurements of the spectra are given in lSergeev et al] 
(|l999h . 

The flux calibration was carried out by assuming the 
narrow emission line fluxes to be constant. We chose the nar- 
row [OIIIJA5007 (for the H/3 region) and [SII]AA6717, 6732 
emission lines (for the Ha region) as the internal flux 
standards. Their absolute fluxes were measured from the 
spectra obtained under photometric conditions by us- 
ing the spectra of the comparison star BS 3082. The 
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Figure 1. Mean Mrk 6 spectrum obtained by combining the 
quasi-simultaneous pairs of spectra from the Ho and H/3 spectral 
regions (SPEC-10 CCD only). 



mean fluxes in these lines are given in ISergeev et al.l 
(1999): F([OIII]A5007) = (6.90±0.11) x KT 13 ergs" 1 cm" 2 ; 
F[SII] AA(6717 + 6731)=(1.61±0.09) x lCT 13 ergs" 1 cm" 2 , 
and F[OI] A6300=(0.597±0.035) x 10 -13 erg s _1 cm" 2 . The 
[SII] lines reside on the far wings of the broad Ha line. Thus, 
to measure their fluxes, we selected the pseudo-continuum 
zones closely spaced around each line. The line fluxes were 
measured by integrating the spectra over the specified wave- 
length intervals and above the continuum (or local pseudo- 
continuum), which was fitted with a straight line in the se- 
lected zones. The mean continuum flux per unit wavelength 
was determined in two windows: at 5162-5186 A (desig- 
nated as _F5170) and 6985-7069 A (designated as F7030). 
T he continuum zones a nd integration limits are the same as 
in lSergeev et al.l (| 19991 ). The line and continuum flux uncer- 
tainties contain errors related to the S/N of the source spec- 
tra, atmospheric dispersion, changes in the position angle of 
the slit, and se eing effects. Evaluati on of these uncertainties 
is considered in lSergeev et al.l (1999). The mean spectrum of 
Mrk 6 produced by combining 23 quasi-simultaneous pairs 
of spectra from the Ha and H/3 regions is shown in FigfT] 
Figure [5] shows the mean and rms spectra of Mrk 6 based 
on our observations in 1992-2008. 



2.2 Optical photometry 

In order to improve the time resolution of our data set 
in the F5170 continuum we added the V-band photome- 
try to our spectral observations. Photometric data came 
from two sources: the UBV observations were made at the 
Crimean Laboratory of the Sternberg Astronomical Insti- 
tute of Moscow University, and the BVRI observations 
were obtained at the CrAO. The UBV observations were 
obtained in the standard Johnson photometric system and 
were carried out from 1986 to 2009 at the 60-cm Zeiss tele- 
scope with a photo-multiplier detector through the aperture 
A=27'.'5. The mean uncertainty of Mrk 6 in t he 1^-band is 
0^022. These data were partially published by iDoroshenkd 
|2003t ). 

In 2001 we started regular observations of Mrk 6 using 
the CrAO 70-cm AZT-8 telescope and the AP7p CCD. The 
CCD field covers 15' x 15'. Photometric fluxes were measured 
within an aperture of 15'.'0. The mean uncertainty of the 
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Figure 2. Mean and rms spectra of Mrk 6 in the H/3 region from 
our observations 1992-2008. 



V-band CCD observations is 0^009. Further details about 
the instrumentation, reductions, and measurements of the 



BVRI photometric data can be found in Doroshcnko et al 



(2005). These data were partially published bv lSergeev et al 

fcooa ). 

2.3 Light curves 

Light curves in Ha and H,9, and the adjacent continuum 
are shown in Fig. [3] for spectral observations corrected for 
seeing. The continuum light curves obtained from the pho- 
tometric IZ-band observations were scaled to the flux den- 
sity measured from the spectroscopic observations. To this 
end, we used the observations made on the same nights al- 
most simultaneously. We have 29 appropriate observational 
nights at the 2.6-m and 70-cm telescopes (spectra plus CCD 
photometry) and 40 appropriate nights at the 2.6-m and 
60-cm telescopes (spectra plus the UBV photoelectric pho- 
tometry). The correlation coefficient between the spectral 
continuum and the l/-band CCD flux for the appropriate 
nights is r=0.984 (n=29 points), and the correlation coeffi- 
cient between the spectral continuum and the U-band flux 
from the UBV observations is r=0.993 (n=40 points). Using 
the regression equations, we converted our U-band photo- 
metric fluxes to the spectral continuum fluxes (F5170). For 
the nights where both spectral and photometric observa- 
tions were available, the continuum fluxes are calculated as 
the weighted average. The fluxes were not corrected for the 
host starlight contamination and Galactic reddening. 

Tables [l] and [2] give the light curves from the 
spectral observations for 1998—2008 (the da ta for 
1991-1997 are available in lSergeev et al.lll999h . The 
H/3 and H7 line fluxes together with the spectral 
continuum fluxes F5170 are shown in Table [Tl The Ha 
line fluxes and the spectral continuum fluxes -F7030 
are shown in Table [2j Table l3l gives combined contin- 
uum fluxes from both the spectral observations and 
from V-band photometric measurements for 1991— 
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Table 1. F5170 spectral continuum, the H/3 and H7 line fluxes. 
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Figure 3. Light curves of Mrk 6 shown for the continuums and 
Ha, H/3, and H7 lines. Narrow lines were not subtracted from the 
fluxes. Bottom panel shows the merged continuum light curve 
used for cross-correlation analysis. Units are 10 — 13 erg cm -2 s _1 
and 10 — 15 erg cm -2 s _1 A -1 for the lines and continua, respec- 
tively. The vertical dash-dotted lines show the boundaries of the 
five time intervals considered in the present paper. 



2008. All the fluxes are seeing-corrected. These light 
curves have been used for the subsequent time-series 
analysis. 

Bottom panel in Figure [3] shows the combined light 
curve from different telescopes. The combined light curve 
shows long time scale continuum variability in Mrk 6 as well 
as more rapid random changes. The flux maxima were ob- 
served in 1995-1996 and in 2007. 

Statistical parameters of the light curves for the lines 
and continuum are listed in Table [4] Column 1 gives the 
spectral features; column 2 gives a number of data points; 
column 3 is the median time interval between the data 
points. The combined continuum light curve is sampled bet- 
ter than the H/3 light curve, the H/3 light curve is sampled 
better than Ha. The mean flux and standard deviation are 
given in columns 4 and 5, and column 6 lists the variance 
F var calculated as the ratio of the rms fluctuation, corrected 
for the effect of measurement errors, to the mean flux. R m ax 
in column 7 is the ratio between the maximum and minimum 
fluxes corrected for the measurement errors. Uncertainties in 
F var and R m ax were computed assuming that a light curve 
is a set of statistically dependent values, i.e., a random pro- 



JD-2,440,000 


F5170 


H/3 H7 


10869.449 


5.932 ±0.057 


3.980 ± 0.079 


10875.449 


6.279 ±0.087 


3.987 ± 0.093 


10876.379 


6.365 ± 0.085 


4.082 ±0.112 


10905.395 


6.451 ±0.057 


4.420 ± 0.080 


10906.273 


6.352 ±0.088 


4.329 ± 0.098 


10924.432 


5.684 ± 0.079 


4.511 ±0.085 


10966.379 


5.008 ± 0.058 


3.476 ± 0.065 


10982.395 


5.570 ±0.101 


3.792 ± 0.096 


10994.387 


5.678 ±0.062 


4.072 ± 0.085 


11014.516 


5.302 ±0.091 


3.971 ± 0.094 


11052.563 


5.519 ±0.081 


3.724 ±0.100 


11074.574 


5.665 ± 0.064 


3.939 ± 0.086 


11076.516 


5.628 ± 0.080 


4.022 ± 0.078 


11085.570 


6.117 ±0.066 


4.102 ± 0.087 


11100.512 


6.344 ± 0.076 


4.283 ± 0.079 


11141.406 


5.155 ±0.087 


4.316 ±0.087 


11218.281 


5.525 ±0.117 


4.022 ±0.133 


11278.395 


6.743 ± 0.074 


4.795 ±0.103 


11281.383 


6.927 ±0.087 


4.687 ± 0.099 


11290.402 


5.878 ±0.187 


4.668 ±0.201 


11310.289 


5.216 ±0.067 


4.594 ± 0.090 


11319.313 


4.817 ±0.063 


4.276 ± 0.086 


11322.305 


5.134 ±0.148 


3.885 ±0.113 


11349.324 


4.884 ± 0.084 


3.615 ±0.078 


11485.504 


5.096 ± 0.072 


4.438 ± 0.089 


11497.434 


5.160 ±0.086 


4.141 ±0.113 


11516.551 


5.103 ±0.073 


4.020 ± 0.078 


11557.418 


4.465 ± 0.077 


3.188 ± 0.098 


11577.406 


4.456 ± 0.066 


3.151 ±0.087 


11587.293 


4.514 ±0.070 


3.119 ±0.079 


11606.383 


4.821 ±0.112 


3.192 ±0.127 


11608.285 


4.540 ±0.089 


3.187 ±0.094 


11615.340 


5.173 ±0.114 


3.081 ±0.108 


11636.270 


5.252 ±0.076 


3.561 ± 0.084 


11707.402 


5.783 ±0.095 


3.687 ±0.092 


11720.336 


5.875 ±0.069 


3.671 ± 0.087 


11780.555 


5.757 ±0.102 


4.052 ±0.123 


11782.555 


5.868 ± 0.080 


3.989 ±0.101 


11791.531 


5.385 ±0.093 


4.017 ±0.126 


11810.531 


5.463 ± 0.075 


3.751 ± 0.097 


11821.582 


5.010 ±0.092 


4.032 ± 0.094 


11823.496 


5.049 ± 0.072 


4.068 ± 0.097 


11838.508 


5.234 ±0.183 


3.741 ±0.181 


11840.539 


5.027 ±0.103 


3.574 ±0.111 


11844.496 


4.679 ±0.072 


3.577 ±0.098 


11847.473 


4.688 ± 0.089 


3.433 ±0.116 


11853.582 


4.628 ±0.067 


3.456 ± 0.089 


11869.496 


5.184 ±0.073 


3.363 ± 0.082 


11878.449 


5.012 ±0.080 


3.450 ±0.102 


11901.316 


5.436 ±0.106 


3.354 ±0.123 


11926.270 


4.391 ±0.084 


3.374 ± 0.080 


11959.227 


4.689 ±0.129 


3.360 ±0.151 


11999.270 


4.517 ±0.112 


3.077 ±0.103 


12030.422 


5.191 ±0.100 


3.149 ±0.105 


12050.465 


4.689 ± 0.087 


3.219 ±0.081 


12104.316 


4.394 ± 0.069 


2.805 ± 0.079 


12174.539 


4.327 ±0.069 


3.007 ± 0.087 


12175.531 


4.153 ±0.084 


3.094 ± 0.087 


12192.551 


4.046 ±0.095 


3.099 ±0.104 


12201.598 


3.845 ±0.105 


2.994 ±0.133 


12223.543 


4.004 ± 0.085 


2.934 ±0.107 


12323.375 


3.899 ± 0.054 


2.424 ± 0.063 


12343.473 


3.573 ±0.107 


2.723 ±0.118 
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Table 1 - continued F5170 spectral continuum, the H/3 and H7 
line fluxes. 



Table 1 - continued F5170 spectral continuum, the H/3 and H7 
line fluxes. 



JD-2,440,000 


F5170 


H/3 


H7 


12381.262 


4.039 ±0.125 


2.356 ±0.130 


- 


12441.359 


5.067 ± 0.078 


3.406 ± 0.082 


- 


12613.543 


5.260 ± 0.085 


3.500 ±0.099 


- 


12619.543 


5.285 ±0.063 


3.457 ± 0.075 


- 


12675.410 


5.032 ±0.106 


2.854 ±0.126 


- 


12914.605 


5.412 ±0.152 


3.171 ±0.156 


- 


12915.531 


5.223 ±0.067 


3.445 ± 0.092 


- 


13089.324 


4.197 ±0.062 


2.682 ± 0.066 


- 


13097.301 


4.849 ±0.205 


2.603 ±0.111 


- 


13319.496 


3.639 ± 0.065 


2.741 ± 0.055 


- 


13356.555 


3.175 ±0.060 


2.219 ±0.068 


- 


13611.563 


6.791 ±0.081 


4.125 ±0.075 


1.342 ±0.057 


13612.553 


6.812 ±0.062 


4.135 ± 0.058 


1.368 ±0.046 


13621.523 


6.877 ±0.082 


4.264 ± 0.079 


1.417 ±0.064 


13641.566 


6.279 ±0.096 


4.197 ±0.078 


1.347 ±0.069 


13647.570 


6.158 ±0.095 


4.157 ±0.075 


1.260 ±0.062 


13648.594 


6.310 ±0.084 


4.156 ±0.077 


1.332 ±0.069 


13649.559 


6.276 ±0.081 


4.082 ± 0.073 


1.277 ±0.063 


13670.555 


6.284 ± 0.083 


4.020 ±0.068 


1.317 ±0.060 


13676.583 


6.309 ± 0.079 


3.914 ±0.065 


1.235 ±0.055 


13698.469 


6.366 ±0.103 


4.038 ± 0.081 


1.255 ±0.067 


13774.453 


5.677 ±0.189 


4.203 ±0.114 


1.290 ±0.104 


13788.344 


5.997 ±0.091 


4.044 ± 0.067 


1.269 ± 0.058 


13830.238 


6.715 ±0.113 


3.949 ± 0.086 


1.206 ± 0.070 


13875.332 


7.178 ±0.086 


4.282 ± 0.086 


1.264 ±0.071 


13994.566 


7.272 ± 0.080 


4.269 ± 0.077 


1.326 ±0.062 


14013.613 


7.496 ± 0.089 


4.407 ± 0.087 


1.317 ±0.066 


14038.555 


7.741 ± 0.089 


4.177 ±0.086 


1.227 ±0.073 


14048.566 


7.806 ±0.136 


4.069 ±0.100 


1.085 ± 0.082 


14064.551 


7.001 ± 0.085 


4.244 ± 0.081 


1.116 ±0.057 


14065.617 


6.963 ± 0.079 


4.139 ±0.075 


1.195 ±0.063 


14076.461 


6.903 ±0.156 


4.044 ± 0.082 


1.194 ±0.076 


14078.551 


6.923 ±0.101 


3.989 ± 0.090 


1.015 ±0.072 


14086.559 


7.182 ±0.097 


3.924 ± 0.080 


1.062 ± 0.065 


14106.531 


7.897 ± 0.243 


4.107 ±0.186 


- 


14111.473 


7.469 ± 0.079 


4.316 ±0.078 


1.241 ±0.064 


14142.594 


7.459 ± 0.097 


4.418 ±0.090 


1.120 ±0.074 


14156.270 


7.366 ±0.132 


4.429 ±0.102 


1.220 ±0.077 


14331.543 


8.416 ±0.086 


5.376 ± 0.095 


1.647 ±0.073 


14369.516 


7.978 ±0.124 


4.843 ± 0.094 


1.556 ±0.081 


14370.539 


7.390 ±0.197 


4.736 ±0.093 


1.335 ±0.074 


14377.582 


7.570 ±0.110 


4.714 ±0.107 


1.410 ±0.088 


14386.457 


7.490 ± 0.098 


4.184 ±0.097 


1.294 ±0.082 


14391.551 


6.816 ±0.097 


4.384 ± 0.076 


1.327 ±0.061 


14392.539 


7.115 ±0.097 


4.410 ±0.081 


1.331 ±0.066 


14393.504 


7.178 ±0.093 


4.517 ±0.092 


1.396 ±0.071 


14423.559 


6.261 ±0.101 


3.715 ± 0.088 


1.149 ±0.084 


14438.438 


5.934 ±0.100 


3.612 ± 0.075 


1.043 ± 0.066 


14480.414 


7.435 ± 0.082 


4.110 ±0.080 


1.166 ±0.067 


14481.398 


7.342 ± 0.086 


4.108 ±0.083 


1.049 ± 0.065 


14482.434 


7.498 ± 0.091 


4.016 ±0.086 


1.179 ±0.075 


14496.367 


7.517 ±0.114 


4.244 ± 0.094 


1.209 ± 0.074 


14497.418 


7.561 ±0.123 


3.955 ± 0.099 


1.155 ±0.082 


14499.328 


7.615 ±0.090 


4.125 ± 0.088 


1.168 ±0.077 


14508.418 


7.071 ±0.106 


4.272 ± 0.088 


1.258 ±0.076 


14537.340 


6.605 ± 0.099 


3.742 ± 0.086 


1.233 ±0.082 


14554.324 


5.749 ±0.125 


3.736 ±0.099 


0.942 ± 0.083 


14574.309 


6.132 ±0.269 


3.303 ±0.226 


- 


14586.309 


6.284 ± 0.086 


3.550 ± 0.074 


1.064 ±0.069 


14587.281 


6.184 ±0.092 


3.576 ± 0.080 


1.003 ±0.067 


14588.266 


6.264 ± 0.098 


3.577 ±0.082 


1.027 ±0.068 


14590.309 


6.278 ± 0.098 


3.543 ± 0.076 


0.968 ± 0.060 



JD-2,440,000 



F5170 



H/3 



Hi 



14591.273 
14617.344 
14620.301 
14622.426 
14687.547 
14779.629 
14804.566 



6.440 ±0.108 
7.076 ± 0.086 
7.053 ± 0.091 
7.141 ±0.087 
6.498 ± 0.095 
7.584 ±0.127 
8.214 ± 0.093 



3.556 ± 0.086 
3.632 ± 0.081 
3.923 ± 0.090 
3.682 ± 0.084 
4.099 ± 0.073 
4.711 ±0.106 
4.650 ± 0.086 



1.020 ±0.077 
0.955 ±0.062 
1.053 ±0.071 
1.048 ± 0.068 
1.091 ±0.057 
1.328 ± 0.070 
1.305 ±0.070 



Units are 10 13 ergs cm 2 s 1 and 10 15 ergs cm 2 s 1 A 



for 



the lines and continuum, respectively. 



cess. The F var values in Table [4] are the lowest limits of 
actual F V ar because the observed fluxes were not corrected 
for starlight contamination. 



3 CROSS-CORRELATION BETWEEN THE 
CONTINUUM AND THE INTEGRAL 
BALMER LINE FLUX VARIATIONS 

As mentioned in the Introduction, estimating the light travel 
time delay between the continuum and emission line flux 
variations is of special relevance for the determination of 
the BLR size, which, i n turn, can be used for black hole 
mass measurements (see lWandel et al.lll999l ; |Peterson et al.l 
2004). This time delay (or lag) is estimated through the 
cross-correlation function (CCF). Koratkarfc Gaskelll ( 19911 ) 
demonstrated that the CCF centroid gives the luminosity- 
weighted radius, in contrast to the CCF peak, which 
is more in fluenced by gas at small radii, according to 
iGaskell fc Sparkd (jl986u . 

The time delays were compute d using the interpolated 
cross-correlation function (ICCF) ilGaskell fc Sparkd Il986l ; 
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IWhite fc Petersonlll994l ; |Petersonll200ll ). We computed both 
the lag related to the CCF peak (r p k) and the CCF cen- 
troid (r cn ). The CCF centroid was adopted to be mea- 
sured above the correlation level at r Js 0.8r maa: . The lag 
uncertainties were computed using the model-independent 
Monte Carlo flux randomization/r andom subset s e lectio n 
(FR/RSS) technique described by iPeterson et~aH (|l998T) . 
The number of realizations was as large as 4000. The uncer- 
tainties were computed from the distribution function for 
T p k and T crl at the 68% confidence, which corresponds to 
±1<T errors for the normal distribution. 

The spectra for the 1992 season were obtained with the 
2" entrance slit and they were discarded from the CCF anal- 
ysis in order to exclude the aperture effects. The results of 
the cross-correlation analysis for the 1993-2008 interval are 
presented in Table [5] The meaning of symbols "s" and "scp" 
in the first and second columns of Table [S] is the same as in 
Table [4] "s" - from spectral observations only and "scp" - 
from combined spectral and photometric observations ("c" 
- CCD, "p" - photoelectric). 

Figure [4] shows the cross-correlation results for the H/3 
and Ha line fluxes with the continuum as well as autocor- 
relation functions of the continuum for the 1993-2008 time 
interval. Table [5] gives the cross-correlation results. When 
the continuum light curve is a combination of the spectral 
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Table 2. F7030 spectral continuum and H« fluxes. 



JD-2,440,000 


F7030 


Ha 


10869.359 


6.858 ±0.194 


31.959 ±0.796 


10876.453 


7.063 ±0.177 


32.400 ± 0.833 


10905.332 


7.403 ±0.176 


33.322 ± 0.722 


11053.559 


6.227 ±0.177 


30.750 ±0.726 


11075.582 


6.408 ±0.171 


30.532 ±0.713 


11220.219 


6.271 ±0.231 


29.815 ± 1.056 


11279.324 


7.851 ±0.210 


37.009 ± 0.962 


11321.301 


5.697 ±0.265 


33.947 ± 1.258 


11517.488 


5.327 ±0.268 


31.821 ±1.134 


11587.355 


5.174 ±0.195 


27.152 ±0.741 


11638.250 


6.202 ± 0.163 


29.285 ± 0.743 


11707.344 


6.400 ±0.165 


30.001 ±0.712 


11721.324 


6.689 ±0.188 


32.688 ± 0.838 


11823.543 


6.115 ±0.221 


32.067 ±0.951 


11847.531 


5.611 ±0.204 


30.121 ±0.891 


11878.512 


6.134 ±0.216 


30.642 ± 0.918 


11902.461 


6.325 ±0.219 


30.126 ±0.899 


12000.254 


4.590 ±0.381 


28.123 ±1.239 


12031.402 


6.159 ±0.277 


27.928 ±1.015 


12105.324 


5.290 ±0.183 


28.077 ±0.731 


12176.520 


4.917 ±0.209 


27.194 ±0.780 


12193.508 


4.927 ±0.263 


28.147 ±1.009 


12225.438 


3.968 ± 0.598 


28.053 ± 1.987 


12323.457 


4.401 ± 0.166 


22.116 ±0.518 


12343.551 


4.813 ±0.394 


25.054 ±1.615 


12613.594 


5.715 ±0.277 


26.821 ± 1.083 


12675.477 


5.779 ±0.193 


26.451 ±0.827 


13090.277 


4.505 ± 0.297 


24.783 ±0.909 


13356.465 


3.857 ±0.159 


20.999 ± 0.656 


13610.559 


8.018 ±0.215 


31.853 ±0.796 


13621.555 


8.487 ±0.240 


32.350 ± 0.891 


13642.594 


7.377 ±0.235 


31.616 ±1.047 


13648.551 


7.688 ± 0.277 


30.374 ± 1.081 


13648.563 


7.474 ±0.218 


30.213 ±0.859 


13789.316 


7.176 ±0.227 


32.545 ± 0.986 


13831.234 


7.418 ±0.186 


30.618 ±0.825 


13994.586 


8.940 ± 0.228 


34.663 ± 0.875 


14066.457 


8.472 ± 0.257 


34.872 ± 1.008 


14370.563 


8.707 ±0.264 


37.536 ±1.117 


14392.563 


8.377 ±0.206 


36.161 ±0.988 


14393.527 


8.125 ±0.256 


35.540 ± 1.017 


14481.387 


8.255 ± 0.243 


34.970 ± 1.044 


14508.445 


8.117 ±0.377 


36.963 ± 1.866 


14617.324 


7.455 ± 0.205 


31.689 ±0.795 


14620.316 


8.052 ± 0.228 


33.558 ± 0.900 


14622.367 


7.785 ± 0.204 


32.227 ±0.890 


14804.543 


10.004 ± 0.242 


38.158 ±0.965 



Table 3. Combined F5170 continuum fluxes from spectral and 
photometric observations. 



Units are 10 13 ergs cm 2 s x and 10 15 ergs cm 2 s x A x for 
the Hq line and continuum, respectively. 



and the photometric data (designated as F5170scp in Ta- 
ble [5]) then the CCF computation has been carried out with 
the rebinning the F5170scp light curve to the times of obser- 
vations of the first time series. This is because the combined 
continuum light curve has much more data points than the 
continuum light curves which consist from the spectral data 
only (designated as F5170s in Tabled . For the F5170s light 
curve, we follow a standard method of the CCF computation 
with rebinnig both time series. 

The variations of the Ha and H/? fluxes are tightly corre- 
lated as well as variations of the continuum fluxes near both 



JD-2,440,000 


F5170 


Julian Date 


F5170 


8630.5781 


4.072 ±0.107 


9783.2969 


7.666 ± 0.098 


8716.4258 


5.003 ±0.143 


9814.2656 


7.775 ±0.124 


8717.3633 


5.083 ±0.129 


9838.3027 


7.282 ± 0.356 


8927.5117 


4.125 ±0.054 


9839.3398 


6.827 ±0.240 


8983.4009 


4.277 ±0.050 


9867.3389 


7.121 ± 0.303 


9001.2852 


4.483 ±0.151 


9871.3545 


6.683 ± 0.252 


9031.2314 


3.991 ±0.150 


9872.3594 


6.775 ± 0.080 


9057.2441 


4.225 ±0.146 


9980.5746 


6.749 ± 0.058 


9059.3291 


4.483 ±0.171 


10008.5636 


7.303 ± 0.052 


9062.3330 


4.331 ±0.161 


10009.5860 


7.380 ± 0.072 


9070.2891 


4.274 ± 0.068 


10010.5900 


7.445 ± 0.083 


9074.3203 


4.378 ±0.059 


10013.5390 


7.236 ± 0.251 


9088.2656 


5.198 ±0.076 


10015.5170 


7.458 ± 0.255 


9089.2578 


5.019 ±0.079 


10024.5360 


7.791 ± 0.289 


9100.2910 


5.021 ±0.186 


10036.5900 


7.798 ±0.119 


9101.2930 


5.116 ±0.174 


10047.5526 


8.101 ±0.057 


9141.3047 


6.684 ±0.119 


10064.5120 


7.593 ± 0.072 


9156.3203 


7.048 ±0.126 


10069.4550 


6.782 ± 0.231 


9250.5977 


6.075 ±0.170 


10092.2260 


7.053 ± 0.241 


9252.5757 


6.082 ± 0.093 


10094.2990 


7.399 ±0.253 


9255.4814 


5.841 ±0.210 


10096.3280 


6.939 ±0.234 


9272.6016 


5.799 ± 0.068 


10102.3720 


7.282 ± 0.255 


9273.4551 


5.769 ±0.211 


10133.2640 


6.606 ±0.229 


9274.5796 


5.568 ± 0.047 


10135.3040 


7.540 ±0.276 


9275.4990 


5.164 ±0.184 


10139.4580 


7.376 ± 0.252 


9311.5742 


5.665 ± 0.077 


10156.3134 


7.602 ± 0.067 


9313.4385 


5.567 ±0.061 


10159.3190 


7.949 ±0.275 


9329.4219 


5.841 ±0.214 


10161.4450 


7.790 ± 0.082 


9331.4531 


6.418 ±0.068 


10201.5550 


6.994 ±0.165 


9332.4521 


6.027 ±0.212 


10202.4060 


7.008 ± 0.092 


9341.3809 


6.313 ±0.219 


10212.2930 


5.893 ± 0.229 


9357.3135 


6.540 ±0.233 


10213.4450 


6.383 ± 0.078 


9359.3467 


6.090 ±0.211 


10218.3130 


5.790 ±0.225 


9362.3525 


6.017 ±0.208 


10222.3060 


6.101 ±0.223 


9364.3281 


6.334 ±0.067 


10225.3320 


5.847 ±0.077 


9365.4180 


6.448 ± 0.060 


10246.3710 


5.474 ± 0.060 


9395.3750 


6.140 ±0.068 


10258.3980 


5.091 ±0.110 


9399.3984 


6.155 ±0.088 


10304.4220 


6.479 ± 0.081 


9450.3438 


6.618 ±0.062 


10361.5310 


6.989 ± 0.068 


9452.2461 


6.651 ±0.100 


10363.4840 


7.053 ± 0.070 


9454.2734 


6.731 ±0.077 


10364.5410 


7.064 ±0.274 


9488.3398 


7.461 ± 0.072 


10372.4804 


7.427 ± 0.040 


9520.3281 


6.829 ±0.076 


10392.5060 


7.411 ±0.271 


9536.3867 


6.575 ±0.064 


10395.4687 


7.013 ± 0.049 


9548.5000 


6.048 ± 0.099 


10396.4610 


7.103 ± 0.063 


9554.5234 


5.900 ±0.111 


10397.4610 


7.274 ± 0.063 


9555.5117 


5.813 ±0.124 


10399.4490 


7.179 ± 0.078 


9566.4570 


5.304 ± 0.068 


10401.4600 


7.018 ± 0.246 


9578.3984 


5.070 ± 0.058 


10403.4010 


7.599 ±0.267 


9599.3672 


4.653 ±0.134 


10404.4920 


7.635 ± 0.258 


9622.2734 


5.361 ±0.094 


10406.5080 


7.446 ± 0.240 


9639.5547 


5.878 ±0.096 


10408.5900 


7.202 ± 0.229 


9653.5635 


6.298 ± 0.074 


10430.3770 


6.573 ± 0.224 


9665.5331 


5.784 ± 0.069 


10434.3200 


6.386 ± 0.055 


9685.3984 


6.646 ±0.147 


10435.3240 


6.586 ± 0.072 


9691.5078 


6.915 ±0.086 


10436.4100 


6.371 ±0.052 


9713.4539 


7.374 ± 0.052 


10461.5560 


5.841 ±0.195 


9716.5488 


7.064 ±0.252 


10482.5040 


6.621 ± 0.088 


9723.4297 


7.226 ±0.081 


10483.4790 


6.683 ± 0.220 


9744.3984 


8.089 ± 0.078 


10484.4480 


6.628 ± 0.209 


9753.3086 


8.467 ± 0.074 


10487.3660 


6.410 ±0.213 


9754.3789 


8.249 ±0.183 


10491.4380 


6.573 ± 0.222 
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Table 3 - continued Combined F5170 continuum fluxes from 
spectral and photometric observations. 



Table 3 - continued Combined F5170 continuum fluxes from 
spectral and photometric observations. 



JD-2,440,000 


F5170 


Julian Date 


F5170 


JD-2,440,000 


F5170 


Julian Date 


F5170 


9771.3086 


7.988 ±0.071 


10495.3670 


6.383 ± 0.057 


11074.5740 


5.665 ± 0.064 


11810.5310 


5.463 ± 0.075 


10509.3810 


6.054 ± 0.062 


11085.5700 


6.117 ±0.066 


11076.4995 


5.650 ±0.077 


11817.5320 


5.986 ± 0.243 


10510.4262 


5.990 ± 0.044 


11088.4720 


5.330 ±0.184 


11818.5550 


5.408 ±0.185 


12349.3890 


3.483 ±0.133 


10511.3540 


5.608 ±0.196 


11100.5120 


6.344 ± 0.076 


11821.5820 


5.010 ±0.092 


12366.2620 


3.890 ± 0.024 


10518.4320 


5.408 ±0.176 


11105.5480 


6.584 ±0.214 


11823.4960 


5.049 ±0.072 


12367.2780 


3.984 ± 0.026 


10519.4260 


5.638 ±0.190 


11110.5780 


5.790 ±0.195 


11838.5080 


5.234 ±0.183 


12368.3450 


4.077 ±0.057 


10521.3830 


5.709 ±0.192 


11111.5760 


5.924 ± 0.207 


11840.5390 


5.027 ±0.103 


12369.2769 


4.188 ±0.022 


10522.3851 


5.517 ±0.067 


11141.4334 


5.127 ±0.081 


11842.4730 


5.388 ± 0.206 


12381.2620 


4.039 ±0.125 


10541.3549 


5.524 ± 0.046 


11163.4850 


5.729 ±0.185 


11843.5660 


5.184 ±0.439 


12385.3530 


4.388 ± 0.033 


10543.2960 


5.059 ±0.188 


11164.3130 


6.207 ±0.220 


11844.4960 


4.679 ±0.072 


12386.2840 


4.481 ± 0.038 


10566.3030 


6.464 ± 0.230 


11176.5790 


5.478 ±0.192 


11847.5150 


4.600 ±0.062 


12387.2960 


4.401 ± 0.034 


10569.3380 


5.517 ±0.196 


11192.3790 


5.810 ±0.201 


11853.5681 


4.767 ± 0.049 


12388.3290 


4.504 ±0.050 


10574.3634 


5.080 ± 0.038 


11197.2630 


5.934 ±0.211 


11866.4100 


5.478 ± 0.233 


12399.3050 


4.898 ±0.169 


10575.3160 


4.840 ±0.054 


11199.3480 


5.021 ±0.174 


11867.3210 


5.155 ±0.197 


12403.3180 


5.097 ±0.190 


10576.3520 


4.868 ± 0.058 


11218.2810 


5.525 ±0.117 


11868.3710 


5.349 ±0.190 


12404.2690 


5.073 ± 0.029 


10580.3820 


5.252 ±0.210 


11261.4920 


7.018 ±0.250 


11869.4960 


5.184 ±0.073 


12405.3055 


5.015 ±0.030 


10597.3852 


5.460 ± 0.083 


11274.3120 


6.951 ± 0.229 


11878.4490 


5.012 ±0.080 


12406.3010 


5.193 ±0.227 


10601.3550 


5.540 ±0.064 


11278.3950 


6.743 ± 0.074 


11879.2630 


5.222 ±0.212 


12407.2730 


4.971 ± 0.034 


10611.3440 


6.296 ±0.099 


11279.2870 


6.650 ±0.237 


11882.4800 


4.656 ±0.170 


12408.2942 


5.104 ±0.025 


10628.3590 


6.778 ±0.081 


11281.3662 


6.943 ±0.081 


11901.3160 


5.436 ±0.106 


12409.3480 


5.078 ± 0.216 


10642.3320 


6.824 ± 0.068 


11290.4020 


5.878 ±0.187 


11902.2690 


5.688 ±0.205 


12410.2870 


5.132 ± 0.028 


10654.3320 


6.808 ±0.153 


11306.3890 


5.709 ± 0.203 


11912.5120 


5.116 ±0.190 


12411.2860 


5.119 ±0.029 


10655.3280 


6.764 ±0.081 


11310.2890 


5.216 ±0.067 


11926.2770 


4.424 ± 0.056 


12417.2820 


5.289 ± 0.052 


10687.5700 


6.033 ±0.116 


11319.3130 


4.817 ±0.063 


11932.4920 


4.795 ±0.184 


12419.3020 


5.272 ± 0.049 


10697.4208 


5.833 ± 0.048 


11322.3050 


5.134 ±0.148 


11959.2270 


4.689 ±0.129 


12421.2880 


5.175 ±0.058 


10699.5270 


5.661 ±0.057 


11346.3850 


5.164 ±0.203 


11999.2700 


4.517 ±0.112 


12440.2860 


5.165 ± 0.050 


10705.5590 


5.587 ±0.188 


11349.3240 


4.884 ±0.084 


12030.4220 


5.191 ±0.100 


12441.3590 


5.067 ± 0.078 


10714.5080 


5.504 ±0.063 


11400.5350 


4.851 ± 0.247 


12050.4650 


4.689 ± 0.087 


12442.2850 


5.123 ± 0.044 


10715.5200 


5.526 ±0.061 


11407.5290 


5.021 ±0.167 


12104.3160 


4.394 ± 0.069 


12459.2830 


5.198 ±0.064 


10728.5390 


5.333 ±0.052 


11409.5370 


4.841 ±0.191 


12139.5420 


4.138 ±0.196 


12476.2920 


5.265 ± 0.063 


10729.5080 


5.399 ±0.056 


11454.5560 


5.658 ± 0.207 


12144.5390 


4.620 ± 0.200 


12484.3370 


4.793 ± 0.064 


10747.5430 


5.962 ±0.076 


11467.4890 


5.310 ±0.181 


12147.5200 


4.060 ±0.155 


12498.4950 


4.379 ± 0.024 


10748.4650 


5.950 ± 0.057 


11485.5040 


5.096 ± 0.072 


12166.5140 


4.629 ±0.167 


12530.5830 


4.424 ±0.031 


10755.5340 


6.196 ±0.211 


11488.6040 


5.136 ±0.193 


12174.5390 


4.327 ±0.069 


12536.5580 


4.640 ±0.025 


10758.5249 


5.924 ±0.052 


11493.3450 


5.438 ± 0.202 


12175.5310 


4.153 ±0.084 


12539.5520 


4.732 ± 0.040 


10759.4380 


5.777 ± 0.076 


11497.4340 


5.160 ±0.086 


12192.5510 


4.046 ±0.095 


12541.5710 


4.653 ± 0.041 


10760.5280 


5.924 ±0.195 


11516.5510 


5.103 ± 0.073 


12199.5620 


4.094 ±0.151 


12557.5940 


4.579 ± 0.027 


10761.5130 


6.185 ±0.206 


11522.5520 


4.610 ±0.162 


12201.5449 


3.888 ±0.054 


12566.4570 


4.422 ± 0.036 


10762.5670 


6.196 ±0.209 


11524.5720 


5.184 ±0.231 


12210.5970 


4.129 ±0.149 


12569.5400 


4.648 ± 0.043 


10777.4300 


5.776 ±0.067 


11525.4560 


5.126 ±0.190 


12223.5055 


4.090 ± 0.060 


12593.6361 


4.954 ±0.035 


10801.2660 


5.955 ±0.203 


11557.3624 


4.581 ± 0.056 


12225.3400 


4.340 ±0.161 


12595.3070 


4.844 ±0.035 


10817.3570 


5.389 ±0.177 


11577.4060 


4.456 ± 0.066 


12231.5060 


4.060 ±0.151 


12596.5600 


4.890 ± 0.043 


10863.3710 


6.562 ±0.221 


11581.3890 


4.547 ±0.153 


12263.4680 


4.094 ± 0.024 


12597.5780 


4.956 ± 0.044 


10866.3790 


6.377 ±0.234 


11586.4480 


4.601 ±0.157 


12265.5230 


4.182 ±0.140 


12605.5350 


5.416 ±0.051 


10867.4230 


6.291 ±0.210 


11587.3295 


4.542 ± 0.053 


12280.5470 


3.783 ± 0.020 


12608.5690 


5.428 ± 0.034 


10868.4180 


6.453 ±0.212 


11588.5210 


4.907 ±0.163 


12281.4510 


3.725 ±0.018 


12609.5610 


5.467 ± 0.030 


10869.4520 


5.996 ± 0.047 


11598.2560 


4.721 ±0.180 


12283.4430 


3.801 ±0.020 


12610.4815 


5.324 ± 0.028 


10873.4040 


6.410 ±0.208 


11603.2890 


4.879 ±0.192 


12298.4070 


3.623 ±0.022 


12612.5340 


5.550 ±0.030 


10874.2920 


6.324 ±0.228 


11605.4510 


5.145 ±0.191 


12301.3840 


3.625 ±0.028 


12613.4293 


5.497 ± 0.030 


10875.4490 


6.279 ± 0.209 


11606.3830 


4.821 ±0.112 


12307.2820 


3.427 ±0.129 


12614.5590 


5.548 ± 0.033 


10876.3961 


6.207 ±0.063 


11608.3236 


4.493 ± 0.056 


12308.3500 


3.628 ±0.019 


12618.3620 


5.532 ± 0.026 


10905.3950 


6.451 ±0.057 


11612.3560 


5.097 ±0.181 


12309.4384 


3.599 ±0.019 


12619.5547 


5.452 ± 0.026 


10906.2730 


6.352 ±0.088 


11615.3400 


5.173 ±0.114 


12310.4000 


3.307 ±0.128 


12620.4070 


5.375 ± 0.035 


10924.4482 


5.708 ± 0.074 


11628.3200 


5.729 ± 0.209 


12313.4260 


3.307 ±0.132 


12621.4710 


5.477 ± 0.034 


10957.3660 


5.136 ±0.200 


11636.2700 


5.252 ± 0.076 


12314.3730 


3.823 ±0.030 


12625.4390 


5.301 ± 0.036 


10966.3790 


5.008 ±0.058 


11661.2870 


5.310 ±0.194 


12316.4110 


3.906 ±0.020 


12634.2880 


4.730 ±0.164 


10982.3950 


5.570 ±0.101 


11707.4020 


5.783 ± 0.095 


12321.4275 


4.018 ±0.023 


12635.4457 


5.088 ±0.155 


10994.3922 


5.671 ±0.053 


11720.3360 


5.875 ± 0.069 


12322.4080 


4.144 ±0.023 


12636.3986 


5.035 ± 0.042 


11014.5160 


5.302 ±0.091 


11780.5550 


5.757 ±0.102 


12323.3560 


4.108 ±0.021 


12665.3620 


4.728 ± 0.026 


11044.3860 


5.408 ± 0.220 


11782.5550 


5.868 ± 0.080 


12324.3990 


4.152 ±0.126 


12672.3350 


4.974 ±0.038 


11050.5250 


5.498 ±0.181 


11788.5360 


5.638 ± 0.243 


12336.2490 


3.507 ±0.110 


12674.4100 


5.041 ±0.031 


11052.5630 


5.519 ±0.081 


11791.5310 


5.385 ± 0.093 


12342.4700 


3.402 ±0.144 


12675.3846 


5.002 ± 0.028 
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Table 3 - continued Combined F5170 continuum fluxes from 
spectral and photometric observations. 



Table 3 - continued Combined F5170 continuum fluxes from 
spectral and photometric observations. 



JD-2,440,000 


F5170 


Julian Date 


F5170 


JD-2,440,000 


F5170 


Julian Date 


F5170 


12343.4263 


3.559 ± 0.046 


12683.3540 


4.635 ± 0.037 


12984.5580 


4.277 ± 0.034 


13384.4410 


3.389 ± 0.020 


12346.3890 


3.499 ±0.132 


12684.3840 


4.703 ± 0.034 


12985.4810 


4.329 ±0.087 


13410.4190 


3.370 ±0.113 


12348.2790 


3.315 ±0.141 


12685.3890 


4.632 ±0.051 


12988.5170 


4.372 ±0.029 


13411.5060 


3.259 ±0.106 


12689.3540 


4.830 ±0.041 


12998.4690 


4.207 ±0.024 


12996.5188 


4.248 ±0.024 


13412.3370 


2.904 ±0.097 


12694.3330 


4.772 ±0.034 


13003.4349 


4.187 ±0.022 


12997.4890 


4.193 ±0.025 


13419.3710 


3.446 ± 0.047 


12696.3900 


5.087 ±0.195 


13006.4560 


4.084 ± 0.056 


13423.3570 


3.454 ± 0.033 


13822.3430 


6.565 ± 0.043 


12697.3620 


4.836 ±0.025 


13007.4330 


4.046 ± 0.034 


13424.4630 


3.439 ± 0.068 


13823.3110 


6.558 ± 0.044 


12698.3500 


4.913 ±0.027 


13015.4930 


3.971 ±0.057 


13425.3360 


3.493 ± 0.041 


13830.2380 


6.715 ±0.113 


12700.3786 


4.836 ±0.027 


13022.4180 


4.391 ± 0.024 


13434.3090 


3.703 ±0.022 


13837.3880 


6.661 ±0.051 


12701.3560 


4.781 ±0.027 


13023.3250 


4.367 ±0.143 


13436.3070 


3.791 ±0.024 


13839.3400 


6.759 ± 0.068 


12703.2680 


4.711 ±0.175 


13058.3160 


4.776 ± 0.026 


13437.2730 


3.483 ±0.114 


13844.3300 


6.879 ± 0.040 


12710.2760 


4.595 ±0.028 


13071.4320 


4.767 ±0.072 


13441.3348 


3.908 ±0.021 


13845.4410 


6.959 ± 0.047 


12716.2540 


4.167 ±0.033 


13073.3260 


4.521 ± 0.033 


13445.2997 


4.160 ±0.024 


13849.2870 


7.053 ± 0.251 


12722.3690 


4.069 ±0.145 


13077.3130 


4.367 ±0.136 


13446.3470 


4.388 ± 0.042 


13850.2917 


7.117 ±0.041 


12723.4030 


3.889 ±0.142 


13083.3410 


4.438 ± 0.030 


13449.3550 


4.369 ±0.034 


13854.2850 


7.053 ± 0.220 


12724.3289 


4.009 ± 0.024 


13084.3190 


4.500 ±0.029 


13459.3238 


4.875 ±0.026 


13875.3320 


7.178 ± 0.086 


12726.2940 


3.771 ±0.138 


13084.3190 


4.500 ±0.029 


13460.2620 


4.978 ±0.028 


13880.3050 


6.746 ± 0.036 


12727.3100 


4.033 ±0.028 


13085.3510 


4.558 ± 0.028 


13461.2693 


5.067 ±0.027 


13881.2850 


6.795 ± 0.041 


12728.3020 


4.123 ±0.024 


13087.2780 


4.425 ± 0.028 


13462.3190 


5.040 ±0.032 


13953.5550 


6.732 ±0.106 


12729.3010 


4.177 ±0.024 


13089.3240 


4.197 ±0.062 


13464.2240 


5.285 ±0.037 


13959.5530 


6.737 ±0.080 


12730.2840 


4.232 ±0.029 


13097.2904 


4.606 ± 0.041 


13465.3125 


5.363 ± 0.031 


13967.5620 


6.846 ± 0.056 


12739.2820 


4.235 ±0.028 


13098.2500 


4.636 ± 0.038 


13471.2810 


5.882 ± 0.206 


13973.5680 


6.731 ± 0.050 


12740.3020 


4.185 ±0.031 


13105.3080 


4.556 ±0.152 


13476.3240 


5.693 ±0.037 


13986.5780 


7.163 ± 0.090 


12742.3080 


4.205 ±0.035 


13111.3030 


4.978 ± 0.027 


13478.2290 


5.844 ±0.111 


13987.4830 


7.138 ± 0.066 


12744.3070 


4.221 ±0.030 


13112.3115 


4.799 ± 0.025 


13487.2960 


5.574 ±0.029 


13989.5540 


7.280 ±0.091 


12745.2760 


4.272 ±0.037 


13113.3120 


4.864 ±0.028 


13493.2990 


5.485 ± 0.029 


13991.5330 


7.168 ± 0.069 


12751.2799 


4.242 ±0.024 


13114.2551 


5.012 ±0.030 


13495.2680 


5.497 ± 0.034 


13994.5660 


7.272 ± 0.080 


12752.2510 


4.309 ±0.024 


13115.2890 


5.108 ±0.029 


13508.3130 


5.755 ±0.037 


13995.5930 


7.363 ± 0.054 


12754.2770 


4.383 ±0.024 


13117.3010 


4.814 ±0.164 


13509.2980 


5.987 ± 0.045 


14010.5960 


7.574 ± 0.066 


12756.2860 


4.324 ±0.024 


13130.2730 


5.063 ± 0.044 


13611.5630 


6.791 ±0.081 


14013.6130 


7.494 ± 0.089 


12757.2770 


4.296 ±0.024 


13133.2930 


5.075 ± 0.026 


13612.5531 


6.801 ±0.083 


14022.5790 


7.319 ±0.058 


12759.2830 


4.647 ±0.165 


13135.2730 


5.004 ± 0.029 


13621.5230 


6.877 ±0.082 


14023.4520 


7.557 ±0.061 


12766.2680 


4.454 ± 0.025 


13148.2890 


4.770 ±0.030 


13641.5660 


6.279 ±0.096 


14038.5550 


7.741 ± 0.089 


12767.2590 


4.544 ± 0.029 


13149.2750 


4.727 ± 0.045 


13644.5290 


5.884 ±0.039 


14044.6030 


8.011 ±0.123 


12770.2590 


4.320 ±0.038 


13153.2910 


4.802 ± 0.032 


13645.6130 


5.900 ±0.076 


14048.5820 


7.691 ± 0.097 


12771.2560 


4.272 ±0.051 


13154.3280 


4.884 ±0.052 


13647.5700 


6.158 ±0.095 


14059.5450 


7.104 ±0.047 


12774.2750 


4.356 ±0.035 


13277.5570 


5.045 ± 0.079 


13648.6094 


5.937 ± 0.048 


14060.6550 


6.988 ±0.147 


12775.2710 


4.353 ±0.043 


13291.5610 


4.765 ± 0.029 


13649.5913 


6.174 ±0.075 


14062.5450 


7.171 ±0.054 


12778.2980 


4.492 ±0.185 


13292.5770 


4.817 ±0.029 


13650.6170 


5.939 ± 0.054 


14064.5904 


7.159 ± 0.049 


12790.3080 


4.933 ± 0.028 


13296.6160 


4.355 ± 0.027 


13651.5770 


5.914 ± 0.040 


14065.4730 


7.052 ± 0.036 


12791.2820 


4.859 ±0.040 


13300.6100 


4.245 ± 0.029 


13653.6120 


5.885 ±0.031 


14067.5428 


6.865 ± 0.045 


12841.5450 


4.824 ±0.104 


13302.6060 


4.161 ±0.029 


13654.5950 


5.972 ±0.038 


14069.6400 


6.920 ±0.056 


12866.5290 


5.115 ±0.034 


13305.6230 


4.073 ± 0.027 


13670.5550 


6.284 ±0.083 


14076.4610 


6.903 ±0.156 


12883.5410 


5.011 ±0.171 


13307.6120 


4.107 ±0.053 


13676.5860 


6.309 ± 0.079 


14078.5510 


6.923 ±0.101 


12889.5360 


5.408 ±0.169 


13308.6250 


4.040 ± 0.046 


13680.5980 


5.868 ±0.033 


14086.5590 


7.182 ±0.097 


12890.5250 


5.376 ±0.032 


13309.6280 


3.986 ± 0.039 


13683.5360 


5.688 ±0.189 


14091.5530 


7.180 ±0.044 


12903.5440 


5.300 ±0.186 


13313.6020 


4.011 ±0.027 


13698.4690 


6.366 ±0.103 


14106.5310 


7.897 ± 0.243 


12906.5980 


5.109 ±0.031 


13314.5960 


4.109 ±0.027 


13702.4780 


5.872 ± 0.203 


14111.5780 


7.450 ± 0.050 


12907.5800 


4.926 ±0.159 


13315.6030 


4.046 ± 0.024 


13708.5590 


6.684 ±0.036 


14116.4500 


7.638 ± 0.059 


12912.5813 


5.089 ± 0.030 


13317.5910 


3.993 ± 0.027 


13724.5730 


7.009 ± 0.049 


14117.5520 


7.552 ± 0.331 


12913.5490 


5.132 ±0.034 


13318.6095 


3.912 ±0.028 


13728.3990 


6.749 ±0.217 


14118.4140 


7.587 ±0.055 


12914.5768 


5.108 ±0.028 


13319.4960 


3.639 ± 0.065 


13733.5160 


6.918 ±0.039 


14119.4440 


7.704 ±0.051 


12915.5268 


5.260 ± 0.049 


13320.4660 


3.581 ±0.116 


13738.4460 


6.940 ± 0.046 


14121.5630 


7.700 ± 0.047 


12945.6380 


4.830 ±0.047 


13323.5880 


3.991 ±0.133 


13739.5320 


6.813 ±0.047 


14123.4870 


7.720 ±0.064 


12947.5950 


4.664 ±0.032 


13331.5011 


3.818 ±0.021 


13744.5540 


6.880 ± 0.047 


14142.5940 


7.459 ± 0.097 


12965.5343 


4.463 ± 0.022 


13355.4852 


3.397 ±0.019 


13747.4950 


6.735 ± 0.078 


14145.3596 


7.674 ± 0.049 


12966.5492 


4.309 ±0.027 


13356.4777 


3.397 ±0.022 


13749.4240 


6.698 ±0.070 


14146.4190 


7.550 ± 0.044 


12967.5156 


4.435 ± 0.026 


13357.4334 


3.360 ±0.018 


13760.3710 


6.082 ±0.030 


14149.4380 


7.735 ± 0.041 


12968.5385 


4.372 ±0.024 


13358.4690 


3.454 ±0.021 


13761.4160 


6.143 ±0.042 


14150.4260 


7.695 ± 0.038 


12973.4510 


4.025 ±0.132 


13365.5680 


3.342 ± 0.037 


13763.3990 


5.901 ±0.035 


14156.2700 


7.366 ±0.132 


12974.5110 


3.982 ±0.131 


13379.5508 


3.316 ±0.019 


13774.4530 


5.677 ±0.189 


14167.3450 


7.289 ± 0.050 


12983.6310 


4.241 ±0.045 


13383.4821 


3.446 ± 0.019 


13787.3870 


6.090 ± 0.203 


14169.2930 


7.252 ± 0.048 
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Table 3 - continued Combined F5170 continuum fluxes from 
spectral and photometric observations. 



Table 4. Light curve statistics. 



JD-2,440,000 


F5170 


Julian Date 


F5170 


13788.3440 


5.997 ±0.091 


14171.3070 


7.250 ± 0.050 


13790.2820 


6.259 ±0.199 


14174.4000 


7.481 ± 0.243 


13799.2870 


6.168 ±0.046 


14180.3510 


7.516 ±0.260 


13807.3650 


6.250 ±0.057 


14181.2360 


7.412 ± 0.042 


13816.4150 


6.349 ±0.041 


14191.2710 


7.835 ± 0.083 


13820.3430 


6.248 ±0.208 


14200.3340 


8.249 ± 0.058 


14201.2850 


7.755 ±0.269 


14523.3480 


6.185 ±0.197 


14204.2350 


8.256 ±0.058 


14530.3709 


6.735 ± 0.064 


14206.3000 


8.379 ± 0.073 


14532.3500 


6.735 ± 0.044 


14213.2460 


7.978 ± 0.049 


14534.3240 


6.662 ± 0.036 


14220.2400 


8.159 ±0.108 


14535.3485 


6.570 ±0.031 


14234.2640 


8.173 ±0.060 


14536.3470 


6.542 ± 0.038 


14281.2840 


8.751 ±0.203 


14537.3175 


6.452 ± 0.051 


14283.2890 


8.724 ±0.117 


14538.3130 


6.436 ± 0.062 


14331.5430 


8.416 ±0.086 


14542.2750 


6.073 ± 0.056 


14337.5760 


8.331 ±0.118 


14554.3231 


6.021 ± 0.030 


14338.5040 


8.265 ±0.051 


14555.2970 


5.994 ± 0.093 


14369.5160 


7.978 ±0.124 


14565.4980 


5.939 ±0.037 


14370.5390 


7.390 ±0.197 


14567.2910 


5.906 ± 0.034 


14371.4990 


7.464 ± 0.073 


14568.2830 


6.029 ± 0.043 


14372.4860 


7.718 ±0.071 


14574.2652 


6.262 ± 0.052 


14376.5510 


7.373 ±0.061 


14582.4170 


6.398 ±0.035 


14377.5820 


7.570 ±0.110 


14585.2780 


6.093 ± 0.046 


14386.4570 


7.490 ± 0.098 


14586.3090 


6.284 ± 0.086 


14391.6032 


7.187 ±0.051 


14587.2927 


6.252 ± 0.034 


14392.5390 


7.115 ±0.097 


14588.2802 


6.184 ±0.054 


14393.5960 


7.037 ± 0.045 


14590.3090 


6.278 ± 0.098 


14423.5722 


6.184 ±0.037 


14591.2730 


6.440 ±0.108 


14425.6040 


6.198 ±0.032 


14596.3030 


6.331 ±0.101 


14426.5370 


6.223 ±0.047 


14600.2710 


6.832 ± 0.055 


14428.4230 


6.298 ±0.058 


14601.3560 


6.695 ± 0.093 


14434.4860 


6.146 ±0.044 


14602.3840 


6.727 ±0.057 


14438.4380 


5.934 ±0.100 


14603.3040 


6.667 ±0.051 


14439.6060 


6.055 ±0.034 


14604.3070 


6.752 ± 0.056 


14443.5395 


6.118 ±0.031 


14617.2919 


6.978 ± 0.059 


14444.5270 


6.180 ±0.033 


14618.2730 


6.876 ± 0.087 


14465.4540 


7.011 ±0.042 


14620.3010 


7.053 ± 0.091 


14467.5530 


7.210 ±0.043 


14622.4260 


7.140 ±0.087 


14472.4660 


7.411 ±0.041 


14628.3120 


6.876 ± 0.046 


14475.4050 


6.727 ±0.227 


14632.3170 


6.960 ± 0.075 


14476.4683 


7.358 ± 0.048 


14643.3460 


7.117 ±0.044 


14477.4660 


7.463 ± 0.058 


14647.3060 


7.129 ± 0.041 


14478.4300 


7.452 ± 0.053 


14687.5470 


6.498 ± 0.095 


14479.4230 


7.456 ± 0.053 


14713.5470 


7.087 ±0.267 


14480.4251 


7.368 ± 0.049 


14718.5400 


7.288 ±0.051 


14481.4803 


7.412 ±0.054 


14720.5630 


7.213 ±0.227 


14482.4074 


7.451 ± 0.053 


14738.5060 


7.030 ±0.237 


14483.4110 


7.452 ± 0.069 


14740.4520 


7.053 ± 0.241 


14484.3690 


7.183 ±0.069 


14742.5900 


7.451 ± 0.051 


14488.4530 


7.111 ±0.120 


14769.5830 


8.158 ±0.265 


14490.4630 


7.314 ±0.062 


14778.5960 


8.054 ±0.050 


14496.4231 


7.715 ± 0.044 


14779.6064 


8.005 ± 0.041 


14497.4230 


7.654 ± 0.054 


14780.5790 


8.052 ± 0.042 


14498.4680 


7.550 ±0.051 


14781.5400 


7.969 ± 0.048 


14499.4361 


7.641 ± 0.045 


14784.3780 


8.013 ±0.128 


14500.4060 


7.739 ± 0.050 


14786.4900 


8.206 ± 0.084 


14502.4057 


7.653 ±0.052 


14796.4990 


8.174 ±0.050 


14503.4437 


7.605 ± 0.045 


14799.6310 


8.263 ± 0.061 


14508.4223 


7.310 ±0.034 


14801.3500 


8.467 ± 0.064 


14512.5400 


6.981 ±0.131 


14802.5980 


8.346 ± 0.052 


14522.2720 


6.496 ±0.211 


14804.5660 


8.214 ±0.093 



Time 


N 


dt-mcd 


Mean 


std b 


Fvar 


Rmax 


series 

(i) 


(2) 


(days) 
(3) 


Flux 6 
(4) 


(5) 


(6) 


(7) 



1992-2008, JD244S630 2454K04 



F5170s 


235 


14 


F5170scp 


742 


3 


H0 


235 


14 


F7030s 


102 


29 


Hq 


102 


20 



6.093 1.101 0.18±0.03 2.65±0.06 

5.780 1.277 0.22±0.05 3.16±0.08 

4.065 0.658 0.16±0.03 2.70±0.09 

6.897 1.261 0.18±0.04 2.43±0.17 

32.06 3.77 0.11±0.02 1.83±0.09 



a The letters in the first column indicate the origin of the contin- 
uum fluxes: "s" - from the spectral observations only and "scp" 
— from combined spectral and photometric observations ("c" 
CCD photometry, "p" — photoelectric photometry). 
b Mean fluxes and standard deviations (std) are in units of 
10 — 13 erg cm -2 s — 1 and 10 — 15 erg cm -2 s _1 A" 1 for the lines 
and the continuum, respectively. 




T (days) 



x (days) 



Figure 4. CCFs of Mrk 6 between the H/3 and continuum F5170 
fluxes (left panel) as well as between the Ha and continuum F7090 
in 1993—2008 (thick lines). ACFs are shown by dashed lines. 



lines (the correlation coefficients are equal to 0.93 and 0.97, 
respectively). The positive lag values for the 'Ha-H/3' and 
'F7030-F5170' light curves mean that the region of effective 
continuum emission at A5170 is probably smaller than that 
at A 7030 and the region of effective H/3 emission is probably 
smaller than that Ha. The probability that the delay for the 
'F7030-F5170' and 'Ha-H/3' is less than zero equals to 0.023 
and 0.036, respectively. 

Table 5. Cross-correlation results for 1993—2008. 



First set Second set 



Tpk 



Continuum fluxes are in units 10 ls ergs cm 2 s 1 A 1 



H/3 


F5170s 


0.84 


21.5+.^ 


2i.oti;S 


H/3 


F5170scp 


0.83 


21.9+Jj 


2i-o±i3 


Hq 


F7030s 


0.88 


3l.6+_l 4 9 


ls.otf.l 


Hq 


F5170scp 


0.85 


26.8+i° 8 6 


22 0+ 6 ' 9 
zz - u _6.9 


Hq 


H/3 


0.93 


o c+8.6 
S.5_ 4 2 


9 c+0.4 


F7030s 


F5170s 


0.97 


7 n+3.7 
'■^-4.4 


5.5±?:S 


F7030s 


F5170scp 


0.96 


8 4+ 2 ' 1 
°' 4 -3.6 


3-2±fJ 


Ho 


F5170sc 


0.53 


33 1+9.2 
•^ -11.6 


27.ll- 1 


Ho 


H/J 


0.82 


°' M -9.1 


-0.6t 2 8 3 3 



Meaning of symbols "s" and "scp" the same as in Tabld4l 
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Table 6. Time intervals for time series analysis in the H/3 region. 



Table 7. Cross-correlation results for the H/3 line for the five 
subintervals. 



Time 


JD2440000+ 




N 


^■^"tned 


Series 




H/3 


cont a 


H/3 


cont 


1 


09250-09872 


38 


51 sp 


14.0 


8.9 


2 


09980-10777 


54 


90 sp 


11.9 


3.0 


3 


10869-11516 


27 


59 sp 


16.0 


6.1 


4 


11557-13356 


17 


266 scp 


19.1 


2.0 


5 


13611-14804 


58 


242 scp 


10.8 


2.1 



a Meaning of symbols " sp" and " scp" is the same as in Tabld4l 



The time interval 1993-2008, that has been used for 
cross-correlation analysis, is very long. We have divided it 
into five subinterval in order to check whether the lag values 
for individual sub-intervals are the same as for the entire 
1993-2008 period, whether they are changed in time, and 
whether they are correlated with the fluxes and with line 
widths. In particular, the effective region of the broad-line 
emission can depend on the incident continuum flux, so the 
lag can depend on the continuum flux. Under a virialized 
motion of the line-emitting gas the expected relation be- 
tween the lag r and the line width V is as follows: V oc r -1 ' 2 . 
To this end, we used only the H/3 spectra because more reli- 
able lag estimate for this line, and we carried out the cross- 
correlation analysis separately for each of the five time inter- 
vals listed in the Table H3 The first and second int e rvals were 
taken the same as in the paper bv lSergeev et al.l |l999j). 

The cross-correlation results for the five time subinter- 
vals are given in Table [7] For the continuum light curve in 
Table [7] that consists of the spectral data only (i.e., when 
the number of data points in both the H/3 and continuum 
light curves is the same) we have used the standard method 
of the CCF computation with the rebinning both time se- 
ries. When the photometric data are added to the contin- 
uum light curve, only the continuum light curve was re- 
binned to the times of observation of the H/3 light curve. 
As for the entire period, the lag uncertainties for the in- 
dividual subinterv als were computed by using the FR/RSS 
method as given bv lPeterson et al.l i|1998l ). The uncertainties 
for each subinterval were computed using 4000 FR/RSS re- 
alizations, and the probability distributions for both r c „ and 
T p k were calculated. Each of these distributions was found to 
be very different from the normal distribution. To obtain 
the unweighted mean lag and its uncertainties we 
have sequentially (one after another) performed the 
convolution of the five individual distributions and 
then we have scaled the r-axis by dividing it by a 
number of subintervals (i.e., by five). After the con- 
volution operation, the final probability distribution 
was found to be almost normal with the expectation 
and standard deviation as given in Table \7\ for the 
mean lag. 

Since there are large gaps in our time series, we de- 
cided to investigate in more detail their effect on the lag 
uncertainties for Mrk 6. To make sure that the uncertain- 
ties in our lag estimates are realistic, we decided to verify 
the effect of the sampling of our time series to the lag de- 
termination and to compare results with those obtained by 
random subset selection (RSS) method. We have generated 
random time series with the same autocorrelation function 



Subset 


N 


points 


Ten 


' pk 


^max 


Number 


H/3 


F5170 








1 


38 


38 


22 7+ 7 - 4 
zz -'-3.0 


20.5+Ji 


0.939 


2 


54 


54 


20.8+3-0 


18.1+.?;° 


0.942 


3 


27 


27 


io 3+4-0 


i8.2+; : i 


0.761 


1 


17 


199° 


26.2+^ 2 g 3 


14.2±J 4 4 6 


0.898 


5 


58 


230 a 


20.2+1° 


26.8« 3 3 8 


0.794 


Mean value: 




21.4±2.0 


19.3±1.9 




1 


38 


51 6 


21.2+*;° 


21 0+ 4 ' 


0.923 


2 


54 


90 b 


20.7t 3 2 ° 4 


22.014; 3 


0.939 


3 


27 


59 b 


20 - 5 ±7.0 


28.5_ 21 9 


0.683 


1 


17 


266 c 


23.9+^ 7 3 


9.5±- 6 


0.881 


5 


58 


242 c 


20A±li 


21-8+1' 


0.803 


Mean value 




21.1±1.9 


20.5±2.2 





a Continuum light curve was combined from the spectral data 

and from the CCD photometry. 

b Continuum light curve was combined from the spectral data 

and from the photoelectric photometry. 

c Continuum light curve was combined from the spectral data 

and from both the CCD and photoelectric photometry. 



as observed ACF. Stationary random process (or time se- 
ries) with a given ACF can be generated from an array of 
independent random values £1 , £2, £3, . . . , £ n . To do this, it is 
necessary to find a matrix Uij, such that the multiplication 
of the matrix U to the vector £ gives a vector of depen- 
dent random values xi, £2, £3, ■ • • , x n , with a given correla- 
tion matrix ptj = ACF(tj — ti), where ti,t2,ta, . . . ,t n are 
times of observations. The matrix U is related to the corre- 
lation matrix p as follows: 

P=UU T , (1) 

where T denotes a transposition operation. We used our own 
algorithm to compute the matrix U from ACF. 

First we have generated 1000 realizations of the con- 
tinuum light curve with a time resolution of 1 day over a 
period of 6452 days, which is longer than the real observed 
time interval (6175 days). 

To simulate the H/3 light curve from the continuum light 
curve we have experimented with the three kinds of transfer 
functions: (1) delta-function 8(t — 20), (2) n-shaped function 
which is a constant for < r < 40, and (3) triangular func- 
tion which is linearly decreased down to a zero value from 
r — to r = 68.3. Here the lag r are in units of days. After 
the convolution with the simulated continuum light curves, 
all the transfer functions give the H/3 light curve with a lag 
of about 20 days. Next the simulated continuum and line 
light curves were rebinned to real moments of observations 
and the cross-correlation functions were computed for each 
realization. The lag peak and centroid were measured. The 
largest uncertainties were obtained when we used the tri- 
angular transfer function. These uncertainties are only due 



Black Hole Mass in Mrk 6 11 



to the sampling of the observation data. Table [8] gives a 
comparison of both methods for uncertainty estimates (i.e., 
the random time series versus RSS method) for the trian- 
gular transfer function. The columns designated as r p k and 
r cn are a lag with ±la uncertainties computed from the 
random time series, while other two columns designated as 
RSS give ±1(7 uncertainties computed from RSS method for 
T p k and r cn , respectively. The last column of the table is 
the mean CCF peak value obtained from the method of the 
random time series. As can be seen from this table, the RSS 
method gives comparable or larger (up to two times) un- 
certainties than the method of the random time series does. 
The random time series method seems to be more direct way 
to estimate lag uncertainties and since the lag uncertainties 
were computed in this paper by the RSS method, they seem 
to be realistic or slightly overestimated. 

There is a contradiction between our lag estimate 
and the preliminary results on Mrk 6 published in a 
conference proceeding by Grier et al. (NASA ADS tag 
2011nlsg.confE..52G). This new campaign had a nightly 
sampling rate and it spanned 125 nights beginning 2010 
August 31 and ending on 2011 January 3. They claimed 
a lag of 8 ± 3 days for the H/3 line in Mrk 6. We have first 
checked the effect of the removing of linear trends from our 
light curves ( as recommended by some authors, e.g., 
Welsh 1999). We have removed linear trends from the H/3 
and continuum light curves for each subinterval, even if there 
are no such trends exist. Then we have recalculated a mean 
lag value and it was to be 19 days, i.e. two days less. Then 
we have generated random time series (continuum and H/3) 
by exactly the same way as described above, but for the 
sampling rate of Grier et al. in order to check whether the 
duration of the monitoring program is important for the lag 
measurements. We found that with the data sampling of 
Grier et al. our lag estimate must be less by one more day, 
and so a total difference between our and their results must 
be three days. The real difference is much larger than three 
days. An unexpected result of our simulation was very large 
lag uncertainties for the 125-days data sampling and for the 
triangular transfer function for the H/3 line (see above). The 
lag uncertainties were found to be as large as ±6 days! How- 
ever, for the 5(t — 20) transfer function (i.e., when the line 
light curve is simply a shifted version of the continuum light 
curve) the lag uncertainties were found to be as small as 
~ 0.1 days. So, for short-term campaigns and for transfer 
functions with a long tail (i.e., when the line light curve is 
not only shifted, but a strongly smoothed version of the con- 
tinuum light curve), the uncertainties in lag estimates can 
be very large and they are connected to the extrapolation 
of the line and continuum fluxes when computing the CCF, 
not to the data sampling. We concluded that we can only ex- 
plain the difference of three days between our and their lag 
measurements. Probably, the rest of the difference is due to 
real changes of lag or due to the measurement uncertainties 
and their underestimation. 

It can be seen from Table [8] that the simulated lag val- 
ues are almost the same for all subinterval. The expected 
lag value can be obtained by convolving the ACF with the 
transfer function and it was found to be: r p k — 20.7 days 
and Ten = 22.4 day in an excellent agreement with the sim- 
ulation results. So, the large gaps in our time series do not 
shift the lag measurements. 



Table 8. Comparison of the lag uncertainty estimates between 
the method of the random time series and the RSS method. 

Interval T p k RSS r cn RSS r s im 

1 20.8±2'i 2.1 22.2±2^ 2.7 0.968 

2 20.8+^ 1.6 22.2+J'g 2.1 0.960 



20.5"* 



6.0 22.4+t? 4.3 0.940 



4 20.4tU 3.4 22.ll*'o 8.5 0.980 

3.3 0.978 



20.6+27 6 - 9 22A - 



Table 9. The broad H/3, and Ha line width measurements for 
1993-2008. 





FWHM 

kms 


^lirte 
-1 


H/3(mcan-spcctrum) 
H/3 (RMS-spectrum) 


6278 ± 378 
5445 ± 468 


2821 ± 13 
2884 ± 97 


Ha (mean-spectrum) 
Ha(rms-spectrum) 


5322 ± 142 
4443 ± 338 


2870 ± 22 
2780 ± 35 



4 LINE WIDTH MEASUREMENTS 

It is well known that the emission-line profile evolution 
cannot be entirely attributed to the reverberation effect 
and that the profile changes usually occur on a time scale 
that is much longer than t he flux- variability time scale (see 
IWanders fc Peterson! Il996i ). To decrease the effect of the 
long-term profile changes on the line width measurements 
and to get sufficient statistics, we measured the H/3 and Ha 
line widths for the five subintervals presented in Table [6] 
The line width is typically characterized by its full width 
at half maximum (FWHM) or by the second moment of 
the line profile, denoted as au ne . To measure FWHM for 
the mean Ha and H/3 line profiles, we removed narrow lines 
from the broad line profiles. This is not required for rms 
profiles. However, the spectra must be optimally aligned in 
wavelength and in spectral resolution in order to reduce the 
narrow line residuals in the rms profiles. It is difficult to 
measure FWHM because both the mean broad and rms H/3 
profiles are double-peaked, and so the scatter in the FWHM 
measurements is much larger than in ou ne which is well 
defined for arbitrary line profiles. The uncertainties in the 
lin e width were obtained using bootstrap method described 
bv lPeterson et all (|2004T ). The H/3 and Ha line width (both 
FWHM and au„ e ) and their uncertainties are listed in Ta- 
ble [9] for 1993-2008. Table \W\ gives the a Une computed sep- 
arately for each considered subinterval and for the H/3 line 
only. 

We examined the relationship between the H/3 lag, H/3 
width, and the continuum flux (see Fig. |3J. No significant 
correlations among above three parameters were found. In 
particular, the virial relationship between the lag and width 
does not contradict to our data, but the correlation coef- 
ficient between them does not differ significantly from the 
zero value. More subintervals and less lag uncertainties are 
required. 
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Table 10. Results for the five subintervals 



Time Mean Flux™ r cn 

series F(5170) F(H/3) days 

1 2 3 4 



(mean) 
5 



km s 

(rms) 



M 6 D „in units 1O 8 M 



(mean) 
7 



(rms) 
8 



6.446 


3.632 


6.472 


3.976 


5.730 


3.587 


4.594 


2.677 


7.027 


3.523 



21.2 
20.7 



20 - 5 -?:o 



+4.0 
-3.2 

+3.0 
2.4 

+5.6 



23.9 



+ 17.0 



20.4 



+4.6 



2813±13 


2836±48 


2804± 6 


2626±37 


2808±14 


2876±46 


2870±13 


3222±39 


2807± 8 


2864±35 



1.77 
1.72 H 



+0.33 
0.27 



1 70+ ' 46 
1 -' U -0.58 



1.80 
1.50 
1.79" 1 



+0.34 
0.28 

+0.22 



2.07 



+ 1.50 



1.69 



+0.38 



2.62 



1.77 1 



+ 1.80 



Average: 



21.1 ±1.9 2812 ±10 c 



2882 ± 100 c 



1.76 ±0.16 1.85 ±0.21 



a Same flux units as in Table 1 for 5170 A continuum and H/3, respectively. 
6 Using Onken et al (2004) calibration, / = 5.5. 

c The line width and uncertainties were computed as weighted average and assuming 
different expectations of the line widths among individual periods of observations. 
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Figure 5. The relation between the continuum flux and the lag 
(.Ten) for H/3 in Mrk 6, as well as between the H/3 width and the 
lag (Ten). 



5 BLACK HOLE MASS OF MRK 6 

Determination of the black hole mass from reverberation 
mapping rests upon the assumption that the gravity of the 
central super-massive black hole dominates over gas motions 
in the BLR. The black hole mass is defined by the virial 
equation 



Me 



f 



ct(AV) 2 
G 



where r is the measured emission-line time delay, c is the 
speed of light, ct represents the BLR size, and AV is the 
BLR velocity dispersion. The dimensionless parameter / is 
the scaling factor, which depends o n the BLR str u cture , 
kinematics and inclination of BLR. IPeterson et al.l (|2004i ) 
argued that r cn for the time delay r, and (Tune, measured 
from the H/3 emission line in the rms spectrum for the emis- 
sion line width AV, provide the most robust estimates of 
the black ho l e mass with the reverberation technique. Later 
ICollin et al.l (|2006r ) confirmed that in most cases for the 
black hole mass estimate the line dispersion ou ne is more 
suitable than the FWHM, and au ne from the rms-spectrum 
is more suitable than the crime from the mean spectrum. 

We adopt an average value of / = 5.5 based on the as- 
sumption that AGNs follow the same M bh — <r* relationship 
as q uiescent galaxies (IQnken et al.ll2004l ). This is consistent 
with lWoo et al.l (|2010r > and allows easy comparison with pre- 



vious results, but this is about a factor of two larger than the 
value of / computed bv lGraham et al.l i|201ll ). The value of 
/ can be more decreased due to the effect of radiation pres- 
sure, as was explored bv lMarconi et al.l (J200 8. 2009). Mar- 
coni et al. suggested that neglecting the effect of radiation 
pressure can lead to underestimation of the true black hole 
mass, especially in objects close to thei r Edd ingto n limit. 
Discu ssion between iMarconi et all l|2008l . 120091 ) and iNetzerl 
(2009) shows that there are many unclear questions in this 
area. Naturally, a corrective term for radiation pressure will 
decrease the /-factor. 

We calculated the black hole mass for Mrk 6 with the 
use of r cn = 21.1 ± 1.9 for the time delay averaged over five 
time intervals and an ne = 2882 ± 100 kms -1 from the rms 
spectra for H/3. With r cn taken in days and (Al/) in kms -1 , 
and taking into account the time dilation correction for the 
value of T cn , the mass is equal to: 



Mbh/Mq = 0.1952 x / x 



r cn (obs) 



x (Avy 



The black hole mass calculated from the H/3 line is (1.85 ± 
0.21) x 10 8 M Q . For the Ha line, the r c „ = 26.8 ± 7.7 days 
and anne = 2780 ± 35 kms -1 and the black hole mass is 
equal to (2.2±0.6) x 10 8 M Q . 

The black hole masses calculated for each of five periods 
of observations are listed in columns 7 and 8 of Table [TD] for 
the anne from the mean and rms spectra. One can see that 
all the estimates of the black hole mass based on the H/3 line 
are the same within the scatter. 



6 THE BLR SIZE-LUMINOSITY AND 
MASS-LUMINOSITY RELATIONSHIPS 

Many characteristics of the Mrk 6 galaxy are typical for ac- 
tive galaxies. In this connection, it is of interest to see the 
localization of this galaxy on the BLR Radius-Luminosity 
and Mass-Luminosity diagrams. These diagrams determine 
a relationship between fundamental characteristics of AGNs. 
To this end, the Mrk 6 luminosity should be known. Up to 
the present, for many AGNs the luminosity in the rest-frame 
Ao = 5100 A has been corrected for host-galaxy starlight 
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Distance from center, r" 



Figure 6. Surface brightness distribution for Mrk 6 in the V- 
ba nd on the b asis o f multi-aperture photoelectric photometry 
by iNeizvestnvl 1119871) . The solid curve corresponds to the model 

fi(r) = a + br + cr 1 ' 4 . 



contribution within the apertu res used in spectral obser- 
vations (see iBentz et al.ll2009al ). These authors used high- 
resolution Hubble Space Telescope (HST) images to measure 
the starlight contribution. This contribution was found to 
be significant, especially for low-luminosity AGNs. 

We tried to get at least a rough estimate of the 
host-galaxy contr ibution using the observations made by 
INeizvestnvl (| 19871 ) at the Special Astrophysical Observatory 
(SAO) in October 1984 with different apertures from A=4'.'3 
to 55". The surface brightness distribution in the host galaxy 
of the Mrk 6 nucleus calculated on the basis of these mea- 
surements is shown in Fig. [6] 

The galaxy contribution in our 3"xll" spec- 
tral window was found to be V 9a i = 15™6 or 
F 9ai =2.08xlO" 15 ergs" 1 cm" 2 A" 1 . The mean flux 
observed in the continuum near Ao=5100 A is 



F 



{gal-\-nuc) 



=6.093x10" 



ergs^cm^A -1 (see Table [4j 
and, thus, the mean flux corrected for the galaxy con- 
tribution is equal to _F n „ c =4.013xl0" 15 ergs _1 cm" 2 A _1 . 
The variability amplitude F var increases from 18% to 27% 
after accounting for the galaxy contribution. The mean 
flux was also corrected for Galactic reddening accord- 
ing to the NASA/I PAC Extragalactic Database (NED) 
ISchlegel et al.l |l998j). The luminosity was found to be 
AL A (5100) = (2.51±0.78) x 10 43 ergs" 1 adopting the galaxy 
distance D=81 Mpc and when the galaxy contribution is 
removed. 

The bolometric luminosity of the Mrk 6 nucleus was 
adopt ed to be Lboi — 9ALa(5100A) according to lKaspi et al.l 
(2000) and it is equal to L bol (nucl)=2.26 xl0 44 ergs _1 . This 
luminosity is far from the Eddington limit (LEdd), which is 
equal to Z/_Ecid=2.16xl0 46 ergs -1 for a black hole mass of 
1.8xl0 8 M Q . In other words, the Eddington ratio for Mrk 6 
is Lboi/LEdd — 0.01. In this case and because there are no 
clear indications of gas outflow from the BLR, the radiation 
pressure has a negligible effect on the reverberation mass 
estimate. 

The position of the Mrk 6 nucleus on the BLR Size- 
Luminosity diagram is shown in Fig. [7] The BLR size and 
the luminos i ty of o ther galaxies in Fig . [7] are taken from 
IBentz et all (1200931 ) and iDennev et all (Bold). The black 
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Figur e 7. Top: H/3 BLR size vs. luminosity at 5100A according to 
IBentz et aO l|2009ah and lDennev et al.l l|2010] ). The luminosity of 
all nuclei was corrected for host galaxy contribution. The solid line 
is the best fit to the relationship \og(RBLn)=— 21.3+0.519 log(L). 
Bottom: the Mass— Luminosity diagra m. The black hole ma ss of 
the majority of AGNs was taken from I Peterson e t al. (2004J), ex- 
cept for Mrk 290, Mrk 817, NGC 3227, N GC 3516, NGC 4051 , 
and NGC 5548, for which new results by IDennev et al.l J2010I ) 
were used. Solid lines show the Eddington limit L^dd an d its 
10% and 1% fractions. The position of Mrk 6 nucleus is indicated 
on both plots. 



NGC 3516, NGC 4 051, 
IDennev et all (|2010T ). 



for which we used new data from 



hole masses in Fig. [7] are taken from iPeterson et aU l|2004h . 
except for the galaxies Mrk 290, Mrk 817, NGC 3227, 



7 VELOCITY-RESOLVED REVERBERATION 
LAGS 

7.1 Entire time interval: 1993-2008 

The question about whether the direction of gas motion 
can be determined from the response of th e line profile t o 
the continuum changes was firstly raised by iFabrikal (|1980f ) - 
Generally speaking, the BLR gas velocity field can be ran- 
dom circular orbits, radial gas outflow or infall, or Keplerian 
motion. Examples demonstrating how the velocity resolved 
responses can be r elated to differe nt ty pes of BLR gas kine - 
matics are given in lPetersonl |200ll ) and lBentz et al.l (|2009tJ) . 
The random circular orbits generate a symmetric lag profile 
with the highest lag observed around zero velocity. The infall 
kinematics produces longer lags in the blue-shifted emission, 
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Figure 8. Two top panels show the H/3 mean and rms profiles 
(left) and the H« mean profile (right) divided into ten bins of 
equal flux separated by vertical dashed lines. The flux units for the 



mean and rms profiles are 10 ergs cm s 



The third 



panels from the top show the corresponding velocity-resolved time 
lag response, where the delays are plotted at the flux centroid of 
each velocity bin. Here the vertical error bars are la uncertainties 
in the lag for each velocity bin denoted by the horizontal error 
bars. The horizontal dotted and dashed lines on the third pan- 
els mark the mean centroid lag and la uncertainty, respectively, 
for the entire emission line. The bottom panels show the peak 
value of the correlation coefficient between the bin flux in the line 
and continuum. Here the horizontal dotted lines mark the corre- 
lation coefficient r max between the entire emission line and the 
continuum, as calculated in Table [5] 



and the outflow gas produces longer lags in the red-shifted 
emission. 

Home et al. (2004) formulated some important observa- 
tional requirements for determining a reliable velocity field 
of the BLR: (1) the time duration of observations should be 
at least three times larger than the longest timescale of re- 
sponse, (2) the mean time between subsequent observations 
should be at least two times less than the BLR light-crossing 
time, and (3) the velocity sampling AV used for the velocity- 
delay maps should be no less than t h e spe ctral resolution of 
the data. According to lHorne et al.l i|2004f) . such conditions 
can allow one to distinguish clearly between alternative kine- 
matic models of the BLR gas motion. 

In order to obtain the velocity-delay map, we measured 
the lag as a function of velocity in several bins across the line 
profile. We divided both the Ha and H/3 lines into ten bins 
of equal flux, and the width of these bins was no less than 
1000-1500 km s _1 . For each bin we calculated light curves 
from the Balmer line fluxes. Then each of these light curves 
was cross-correlated with the continuum light curve follow- 
ing the same procedure as described in Section [3] Figure [8] 
shows hydrogen line profiles (mean and rms) subdivided into 
bins (two upper panels). The two middle panels demonstrate 
the lag measurements for each of the bins. The vertical error 
bars show la uncertainties for the time lag, and the horizon- 
tal bars represent the bin width. The horizontal solid and 
dashed lines in the two middle panels show the mean BLR 



lag and associated errors as listed in Table [5] The bottom 
panels show the peak correlation coefficient between the bin 
flux in the line and continuum, r m ax. Figure [8] shows that 

(1) The mean and rms profiles of H/3 and Hq are not 
symmetric with respect to zero velocity. The centroid of the 
mean and rms profiles is shifted to the short-wave part of 
the line. The variable parts of H/3 and Hq have two well- 
defined peaks, one of them is almost central (between 5 and 
6 bins) and another is blue-shifted. In addition, there is a 
weaker peak in the red part of the line profile. 

(2) The time delay between the higher velocity gas in the 
BLR and the continuum is shorter than the delay between 
the low velocity gas and the continuum. Such a behaviour 
is typical for virialized gas motions. 

(3) The lag in the blue wing of the H,9 line is greater 
than the lag in the red side of this line. The Ha velocity- 
resolved lags shows the same tendency. This is consistent 
with expectations from the infall model of gas motion. Thus, 
it is possible we have virialized motion combined with infall 
signatures. 

(4) The correlation coefficient of different segments of the 
lines is different. The bin corresponding to a radial veloc- 
ity of V r ~ +1500 fcms -1 shows poor correlation with the 
continuu m variation, especial ly in H/3. This fact was earlier 
noted bv lSereeev et all l|l999l ). 



7.2 Velocity delay maps in the five time intervals 

We have computed the velocity-resolved time delays for the 
five time intervals given in Table [6] For each subset we made 
the velocity-dependent cross-correlation analysis for the H/3 
line profile bins as described in the previous section. The H/3 
line was selected because its sampling is better. In Figure [9] 
the mean and rms H/3 profiles, the velocity-resolved time lag 
response, and the velocity-dependent peak correlation coef- 
ficient are shown for the five time intervals. Upon inspection 
of Figure [9] it becomes clear that the mean and rms profiles 
are different among the five periods. The relative intensity in 
the blue peak and in the central peak changes very strongly: 
during the first interval the blue peak is higher then the red 
one. The opposite situation is seen in the fifth interval. The 
flux in continuum as well the flux in the H/3 line systemati- 
cally decreased from the second to the fourth time interval, 
as is seen in Figure[3]and Figure[7] The H/3 rms profile shows 
two peaks in the first, second and third period, the flat top 
in the fourth period, and in the fifth period we see three 
peaks. 

Figure [8] demonstrates that the high- velocity gas in the 
wings exhibits a shorter lag than the low-velocity gas, sup- 
porting the virial nature of gas motion in BLR: the gas kine- 
matics that is dominated by the central massive object. 

However, the lag is slightly larger in the blue wing than 
in the red wing for all subintervals. This is a signature of 
the infall gas motion. In the fifth period (2005-2008) the 
velocity delay map is more symmetric, but the seventh bin 
shows very small lag, as well as in the previous time in- 
terval. For 2005-2008 there is a poor correlation with the 
continuum for bins 7-10. A more detailed examination of 
the bin light curves for 2005-2008 (Figure \W§ revealed that 
there is a trend for the H/3 flux in bins 7-10, which is al- 
most absent in the bins 1-6. Following the advice of our 



Black Hole Mass in Mrk 6 15 



1.09.20-1 995.06.03 



1995.09.20-1997.11.24 



1998.02.24-1999.12.04 2000.01.13-2002.04.17 



2005.08.29-2008.12.04 



2500 5000 7500 




A 



+ 




^%, 



fcf 




Tf%f 



\ 



HRir 




1.0 
0.8 
0.6 
0.4 



sJr\ 




-2500 2500 5000 -5000 -2500 2500 5000 -5000 -2500 2500 5000 -5000 -2500 2500 5000 -5000 -2500 2500 5000 7500 

Radial velocity (km/s) 
Figure 9. The same as in Figure [8] but for the H/3 line for the five time intervals. 
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Figure 10. Light curves in the continuum and H/3 in 2005-2008 
over the velocity range from —6200 to +4500 km s —1 (top panels) 
as well as the light curves of the broad H/3 line in bin 5 in the 
velocity range from —980 to +25 kms -1 , and in the bin 7 in the 
velocity range from 1000 to 2000 kms" 1 (bottom panels). The 
continuum fluxes are in units of 10 — 15 ergs - 1 cm~ 2 A~ 1 , and 
the broad emission line fluxes are in units of 10 -15 cn-<= — x .-m -2 
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reviewer we removed the trend from the H/3 light curves for 
bins 7-10. No more significant trends were found for other 
time intervals. In Figure [5] the detrended lags and correla- 
tion coefficients are shown by open circles. After detrending 
procedure, the lag-velocity dependence became more simi- 
lar to the lag-velocity dependence for the first period, for 
which the difference in lag between the blue and red wings 
is largest. 

So, it is most likely that the BLR kinematic in Mrk 6 



is a combination of the Keplerian gas motion and infall gas 
motion. 



8 CONCLUSION 

We have reported our new results on the Mrk 6 nucleus from 
1998-2008 observa tional data together with the previous re- 
sults published bv lSergeev et al.l II 19991 ). We found that 

(1) The flux of the Mrk 6 nucleus in 1992-2008 varied 
significantly in the continuum as well as in the Ha and H/3 
broad emission lines. The relative amplitude of the contin- 
uum flux variability is larger than in the hydrogen lines, and 
it is greater in the H/3 line than in the Ha (see Table [4|. It is 
typical for the most of Seyfert galaxie s. This agrees with the 
predictions of lKorista fc Goadl (J200J) based on new photo- 
ionization calculations of the BLR-like gas. 

(2) We found the average time delay between the total 
R/3 flux and the continuum flux at 5170 A to be 21.1 ± 1.8 
days, and the time delay does not vary significant among 
individual time intervals. It seems that the size of the H/3 
emission region remains approximately the same over long 
time periods. The Ha flux responds to the changes in the 
F5170 continuum with a lag of 26.8l4° g 6 days. 

(3) When the continuum flux varies, the photo-ionization 
models predict the existence of the relation between the BLR 
size and the luminosity. However, because the large uncer- 
tainties in the lag for individual time intervals we are unable 
to find such a relation. For the same reason, it is unable to 
obtain a dependence between the lag and line width, and it 
is impossible to check whether this dependence is consistent 
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with the dependence V oc r -1 ' 2 expected for the gravita- 
tionally dominated motion. 

(4) The H/3 line width is larger than that of Ha. This 
is naturally explained by photo- ionization calculations (e.g., 
iKorista fe Goadl 120041 ): the effective emission region of H/3 
is smaller and closer to the ionizing source than the effective 
region of Ha, and the gas velocities in H/3 are higher. 

(5) By examining the velocity-resolved lags for the broad 
H/3 and Ha lines, we found that the lag in the high-velocity 
wings are shorter than in the line core. This indicates virial 
motions of gas in the BLR. However, the lag is slightly larger 
in the blue wing than in the red wing for the entire time in- 
terval as well as for the individual periods considered in the 
present paper. This is a signature of the infall gas motion. 
Probably the BLR kinematic in the Mrk 6 nucleus is a com- 
bination of the Keplerian gas motion and infall gas motion. 

(6) Some profile segments often s how poor cor relation 
with the continuum flux. According to Gaskcll (2010) this ef- 
fect can arise because off-axis sources of ionizing continuum 
flux can appear, which might not make a detectable con- 
tribution to the total continuum flux variability, but they 
will have an influence on the line only over a narrow range 
of radial velocity in the BLR. If these local off-axis events 
will vary out of phase with the variability of the dominant 
source, the result will be to give a weak correlation between 
the continuum flux and the line flux in the narrow range of 
radial velocity. 

(7) We determined the black hole mass from the lag 
and line width measurements of the H/3 and Ha lines. 
The mass was found to be Mbh - (1.8 ± 0.2) x 10 s M Q 
for the H/3 line and slightly greater and less reliable from 
the Ha line. Under such a mass and the luminosity of 
AL A (5100) = (2.51±0.38)xl0 43 ergs -1 , the Mrk 6 nucleus 
is located on the upper edge of the Mass-Luminosity dia- 
gram that corresponds to the Eddi ngton ratio of about 0.01. 
This confirms the assumption ('e.g. JSergeev et al.ll201lr ) that 
there is anticorrelation between broad-line widths and Ed- 
dington luminosity ratio L^ ijLEdd- The Mrk 6 position on 
the BLR Size-Luminosity diag ram does not contra dict the 
fit Rblr oc L - 5 determined bv lBentz et al.1 (|2009aT ). 
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